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We consider quantum excitation energy transport (EET) in a network of two-state nodes in the 
Markovian approximation by employing the Lindblad formulation. We find that EET from an 
initial site, where the excitation is inserted to the sink, is generally inefficient due to the inhibition 
of transport by localization of the excitation wave packet in a symmetric, fully-connected network. 

We demonstrate that the EET efficiency can be significantly increased up to 100% by perturbing 
hopping transport between the initial node and the one connected directly to the sink, while the 
rate of energy transport is highest at a finite value of the hopping parameter. We also show that 
prohibiting hopping between the other nodes which are not directly linked to the sink does not 
improve the efficiency. We show that external dephasing noise in the network plays a constructive 
role for EET in the presence of localization in the network, while in the absence of localization it 
reduces the efficiency of EET. We also consider the influence of off-diagonal disorder in the hopping 
parameters of the network. 

PACS numbers: 03.65.Yz, 03.67.-a 


I. INTRODUCTION 


There are many systems in nature and society that 
can be modelled by classical complex networks and ex¬ 
plained by statistical physicsRecently, the quantum 
world has opened \m new perspectives in the field of 
complex networks.l^^For example, energy, charge, or in¬ 
formation transfer are important phenomena in physi¬ 
cal and biological systems taking place at scales rang¬ 
ing from atoms to large macro-molecular structures, and 
the idea has been put forward that quantum mechan¬ 
ics might have a positive effect on the efficiency of en¬ 
ergy or charge tran^ort in such systems. Charge trans¬ 
port through DNA^ and energy transfer in photosyn¬ 
thetic structure^^K^ are good examples in this context. 
In fact, the most important effect of quantum mechanics 
in biological systems to date has been seen in the Fenna- 
Matthews-Olson (FMO) complexeJ^, observed by exper¬ 
imentally via ultrafast spectroscopy^, where there is an 
ultrahigh efficient excitati on en ergy transport (EET) in 
light-harvesting complexes ililMHI This can be mo deled 
by quantum walks,II^M^ or the Lindblad formalism.E^E^I 
Simulated artificial complex networks based on biologi¬ 
cal systems have the potential to be used in future quan¬ 
tum informational and computational technologies such 
as secure information transfer between two particles in 
an entangled quantum-cry ptography^l teleportation in 
a quantum communication protocol,l2^H^ artificial pho¬ 
tosynthesis systems (or solar cells) to save energy, and 
quantum neural networkJ^^^^ for advanced technologies 
based on artificial intelligence. 

The in fluenc e of edge deletion in networks for quan¬ 
tum stat (j^^ l ^^ and energy transfeJ^ has been studied; 


however, the influence of edge perturbation and deletion 
on the destructive interference of transition amplitudes 
in a networlP^ remains unexplored. In this paper we 
consider this problem in a fully connected network and 
show how the energy can be transferred in the network 
via deletion of a single edge. We analyze quantum en¬ 
ergy transport in a simple model of a complete network 
of two-state nodes. In particular, we demonstrate how 
energy transfer in such a network, which is inhibited by 
localization due to destructive interference, can be en¬ 
hanced by breaking the symmetry between the nodes, 
or introducing disorder in the coupling between the net¬ 
work nodes. Our work is relevant from the point of view 
of building an artificial network for high efficient energy 
transfer to simulate systems such as FMO complexes in 
photosynthetic structures. 


The paper is organized as follows. In Sec. [TTj we de¬ 
scribe the model Hamiltonian of the network and intro¬ 
duce the method of calculating excitation energy trans¬ 
port, employing the Lindblad formulation. In Sec. we 
investigate EET in a fully connected network and discuss 
the negative effect of localization in this case. In Sec. |IV| 
we study the influence of single edge deletions in the net¬ 
work on the efficiency of EET. We then investigate the 
influence of dephasing on EET in Sec. |V| by considering 
edge-deletions. In Sec. |Vl|we examine how the rate of en¬ 
ergy transport can be optimized in the network. Finally, 
in Sec. |VII| we discuss the role of off-diagonal disorder 
in the coupling constants between the network nodes in 
enhancing the efficiency. Section VIII presents our sum¬ 
mary and conclusions. 
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II. QUANTUM EXCITATION ENERGY 
TRANSPORT IN NETWORKS 

Consider a graph G (i.e. a complex network) as a pair 
of sets G = (z^, 5 ), where z/ is a set of vertices or nodes 
of the graph and 5 is a set of edges or links Vij {Vij G e) 
connecting the vertices i and j {i^j G n). An undirected 
graph is completely defined by its adjacency matrix A, 
defined as: 


r 0 , for Vij ^ e; 

I 1; for Vij ^ £• 


( 1 ) 


A complete N-graph (i.e. a fully connected network) is a 
graph with ( 2 ) edges where Vij = 1 for the all pairs of 
nodeJ^. 

Consider a network, consisting of N nodes, in which 
each node is a two-state object such as a molecule or 
a qubit, with a ground and an excited state. We as¬ 
sume that the nodes interact with each other through 
direct hopping. When hopping between two molecules is 
allowed, a link is drawn between them. Whenever an ex¬ 
citation is inserted at one node, it can then be transfered 
throughout the network by hopping due to the interac¬ 
tion between the linked nodes. Here we consider quantum 
energy transport in a network that can be modelled by 
the following tight-binding HamiltoniaiP^ 


N N 

H = y]fiw„|n)(n| + E J nm {\m){n\ + \n){m\), (2) 

n=l n^m 


where |n) is the n-th site in which the excitation exists, 
hwn denotes the excitation energy at site n, and Jnm is 
the hopping integral between the two sites n and m. We 
set Jnm = 1 when the nodes n and m are connected and 
Jnm = 0 if they are disconnected, and we choose h = 1 
so that all the energies are in units of ujn • Possible decay 
events of the exciton to the ground state are neglected 
here. 

To study dissipationless quantum excitation energy 
transport (EET) in the network, we use the evolution 
of the m aster equation in the Markovian approximation 
as followP^^ 


P — P\ “1“ -^sinkP: (3) 

in which p is the density matrix. EginkP is the sink term 
that expresses the irreversible transfer of energy from a 
given node of the network into a sink, and it is defined 
as: 

isinkP = r[2cr+„kC7f- pcr+C7T„k - {crf+<Ts7nk > f'}] > 

(4) 

where the curly brackets denote the anticommutator, and 
r is the absorption rate of the sink which is set to 0.5 
in all of the calculations in this paper. The quantities 
(7^ (cr^) are the creation (annihilation) operators at the 
site connecting site to the sink, and cr^^k (^sink) 
creation (annihilation) operators at the sink. 


Once an excitation is initially injected to the initial 
i-th site, we have p(0) = |'z)(i|. In order to measure 
EET from the initial site to the sink, we integrate the 
master equation ® and calculate the population of the 
sink (sink|p(t)|sin^ at time t. The population of each 
node is {n\p{t)\n). Another quantity of interest is the 
system efficiency which is defined as the long term sink 
population: 


rjoo = lim(_>oo(sink|p(f)|sink), (5) 

and determines the fraction of excitation energy trans¬ 
ferred into the sink in the long time limit. 

Our calculations are done using the python package 
QUTIpl^for numerical integration of the Lindblad mas¬ 
ter equation (§, and all energies, time scales, and rates 
are expressed in the units of on-site exception energies 
ujn, and since we assume that the network iconsists of 
the identical units, we set uJn = I for all n. 


III. FULLY CONNECTED NETWORK 

We first consider EET in a fully connected network. In 
Eig. Eig.j^a) we show a fully connected six-site network 
(ECN). The sink is connected to the node 6 and the dis¬ 
sipationless excitation is injected to the node 1. Eigure 
Bb) shows the time dependence of the population of the 
nodes and the sink, calculated using numerical integra¬ 
tion of Eq. §. It can be seen that the system efficiency 
(sink population) tends to ^ 0.2, while most of injected 
energy (^ 80%) remains inside the network, mostly on 
the first node (~ 64% of energy remains in node 1 and 
~ 16% is shared equally among the nodes 2,3,4 and 5). 
This result is a manifestation of the localization of sin¬ 
gle particle states within a fully connected network as 
already pointed out in Ref. [161 The reason for localiza¬ 
tion is the existence of destructive interference of tran¬ 
sition amplitudes inside the network. To explicitly show 
this, we expand the initial state |1) in terms of the or¬ 
thonormal eigenstates of the tight-binding Hamiltonian 
(i), resulting in: 

|1) = -0.71 (-0.71,0.71,0,0,0,0) 

-0.41 (-0.41, -0.41,0.82,0,0,0) 

-0.29 (-0.29, -0.29, -0.29,0.87,0,0) 

-0.22 (-0.22, -0.22, -0.22, -0.22,0.89,0) 

-0.18 (-0.18, -0.18, -0.18, -0.18, -0.18,0.91) 
+0.41 (0.41,0.41,0.41,0.41,0.41,0.41). (6) 

Eq. clearly shows that in the four out of six terms in 
the expansion of the initial state, the sixth node (which 
connects the network to the sink) has no contribution, 
and hence most of the injected energy cannot reach this 
site and be transferred to the the sink. Eor visualiza¬ 
tion of the energy localization, the long time limit of the 
density matrix of the network is illustrated in Eig. [^c). 
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FIG. 1: (a) Schematic model of a fully connected six-node network (FCN), where site 6 is connected irreversibly to the sink, 

(b) Population at nodes 1, 2, 3,4, 5 and 6 in the network for F = 0.5. Note that the variation of the populations for the nodes 
2, 3, 4, 5 is exactly the same, as expected from symmetry, (c) Density plot of the stationary density matrix p(oo) of the network. 


where the block diagonalization of the stationary den¬ 
sity matrix is an indication of energy localization in the 
network. 

To investigate the effect of symmetry on the energy 
localization, we calculate the system efficiency when the 
hopping integral between the nodes 1 and 6 varies, while 
the rest of hopping integrals Jnm {Tn^n ^ 1,6) remain 
equal to unity. The result is displayed in Fig. and it 
shows that the system efficiency (r^oo) is highly sensitive 
to the value of Jiq. For Jiq = 1, rj^ is minimized, but 
when the hopping integral between the two nodes 1 and 
6 slightly deviates from the other ones, the system ef¬ 
ficiency rapidly rises to unity. This result shows that 
the state localization inside a fully connected network is 
highly sensitive to the symmetry of the Hamiltonian and 
any asymmetry due to Jiq destroys the destructive inter¬ 
ference loops and hence localization in the network. The 
important conclusion here is that such a symmerry break¬ 
ing dramatically increases the efficiency of EET. Eigure 

also shows that introducing asymmetry into the other 
links slightly increases the system efficiency, but does not 
promote the system to become a perfect transmitter. 


IV. INFLUENCE OF EDGE DELETION ON 
SYSTEM EFFICIENCY 


In the previous section we demonstrated that even a 
small asymmetry in the hopping rates is able to destroy 
the state localization in a fully connected network, and 
hence it increases the energy transfer to the sink. Next 
we proceed to investigate the influence of a defect on en¬ 
ergy transport. To this end, we consider cutting one link 
between any two nodes. Deletion of an edge between 
the nodes means blocking the hopping of the excitation 
between them. Figure illustrates the density matrices 
of the networks with deletion of the links between the 
nodes 1,2 (Fig.|^a)), 2,3 (Fig.|^b)) and 1,6 (Fig.j^c)). 
It can be seen from this figure that only in the case of 
deletion of the link between the nodes 1 and 6, the local¬ 
ization vanishes and the injected energy gets completely 



FIG. 2: System efficiency as a function of the value of the the 
hopping integral between the different nodes I and k {Jik)- It 
can be seen that the hopping rate Jie connecting the insertion 
point and the sink has the strongest influence on the efficiency. 


transferred into the sink. Preventing hopping between 
the nodes 1, 2 and 2, 3 results in block diagonalization of 
the steady density matrix and hence localization inside 
the network is preserved. This can be verified by ex¬ 
pansion of the initial state in terms of the eigenstates of 
the Hamiltonian in Eq. (§, which results in the following 
expression; 

|1) = -0.707 (-0.71,0.71,0,0,0,0) 

+0.365 (0.36,0.36,0.43,0.43,0.43,0.43) 

-0.606 (-0.61, -0.61,0.26,0.26,0.26,0.26) ,(7) 

when the link between nodes 1 and 2 is deleted, and 

|1) = -0.707 (-0.71,0,0,0.71,0,0) 

-0.408 (-0.41,0,0, -0.41,0.82,0) 

-0.289 (-0.29,0,0, -0.29, -0.29,0.87) 

+0.428 (0.43,0.36,0.36,0.43,0.43,0.43) 

+0.258 (0.26, -0.61, -0.61,0.26,0.26,0.26) (8) 
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FIG. 3: (Color online) Density plots showing the behaviour 
of the elements of the steady state density matrix, p(oo), in 
the case of deletion of edges between the nodes (a) 1,2; (b) 
2,3 and (c) 1,6. The dotted lines in the left panels denote 
the edges deleted. As can be seen in case (c), localization is 
absent in the network and therefore the EET efficiency to the 
sink becomes unity. 


framework of the Lindblad operators as follow^^^^ 

N 

Ldephp = E]7deph[2CT+cr“/9(7+(7“ - {o-+(7“,p}] (10) 

n=l 

where 7deph is the dephasing rate coefficient. 

Adding Eq. ( p^ to the master equation (§ and in¬ 
tegrating this equation as before gives energy transfer in 
the presence of noise. Figure [^represents the time depen¬ 
dence of the sink population in the absence and presence 
of noise with different values of the dephacing rate co¬ 
efficient. The results show that in cases where the wave 
packet is partially localized in the network, such as in the 
fully connected (Fig.j^a)), 2 — 4 edge deleted (Fig.j^c)), 
and 3 — 6 edge deleted (Fig. |^d)) networks, the sys¬ 
tem efficiency monotonically increases (both in magni¬ 
tude and speed) as dephasing is increased. Therefore, in 
these networks the noise has a constructive influence on 
energy transport by reducing quantum coherence. How¬ 
ever, when the hopping between the initial node and the 
one connecting to the sink (1 — 6) is blocked, the wave 
packet localization is already destroyed and in this case 
the noise slows down the speed of energy transfer to the 
sink as illustrated in Fig. |^b). 


VI. SATURATION TIME 


when the the link between nodes 2 and 3 is deleted. Fi¬ 
nally, 

|1) = -0.707 (-0.71,0,0,0,0,0.71) 

+0.365 (0.36,0.43,0.43,0.43,0.43,0.36) 

0.606 (0.61, -0.26, -0.26, -0.26, -0.26,0.61)(9) 

in the case of cutting the link between nodes 1 and 6. 

Equations 0 and 0 explicitly show the vanishing 
contribution of the node 6 to some of the terms in the ex¬ 
pansion, while Eq. 0 indicates that node 6 contributes 
to all the components of the initial wave packet and in 
this case all the energy can be transmitted to the sink 
through this node. We also checked that edge deletion 
between all the other nodes preserves localization in the 
resulting network. Figure represents the block diago¬ 
nalized structure of the stationary density matrices and 
the existence of localized states for some of these cases. 


V. INFLUENCE OF DEPHASING 

Next we investigate the influence of a noisy environ¬ 
ment on quantum energy transport in netwo rks. In most 
practical cases the networks are open ISMI interact 
with the environment, which is expected to reduce quan¬ 
tum coherence and constructive interference. Dephasing 
effects can be incorporated in the master equation in the 


The results presented here show that the link which 
connects the initial node to the final node plays a fun¬ 
damental role in EET, and as Fig. shows, when Jiq is 
changed from 0.0 to 0.945 or from 1.055 to 2.0, complete 
energy transfer is obtained. An obvious question then 
concerns finding a value of which is optimal in the 
sense that the saturation time Tq is minimized. We de¬ 
fine Ts to be the time after which the injected energy is 
transferred to the sink in the absence of dephasing. In 
Fig. [^we show numerical results for Tq when varying the 
hopping rate Jiq. The nontrivial result here is that ex¬ 
cept for the peak caused by localization, the saturation 
time decreases with increasing Jiq attaining a minimum 
which is numerically determined to be J^q ^3.04. 

VII. INFLUENCE OF OFF-DIAGONAL 
DISORDER 

In addition to edge deletion considered in the previous 
sections, it’s also possible to introduce disorder explicitly 
into the hopping integral parameters between the net¬ 
work nodes. In such networks it is expected that the 
efficiency is strongly affected. We have studied this issue 
by introducing disorder to the couplings by changing J^m 
to Jnm(l + ^): whcrc ^ is a dimensionless random num¬ 
ber with a uniform distribution in the interval [—x,x]. 
In Fig. 7 we show numerical results for the efficiency 
parameter 7]^^ as a function of x for two fully connected 
networks with N = 6 and N = 7 nodes. As expected. 
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FIG. 4: (Color online) Density plots showing the behavior of the elements of the steady state density matrix, p(oo), in the case 
of deletion of edges between nodes (a) 2,4, (b) 2,6, (c) 3,4, (d) 3,6, (e) 4,5, and (f) 5,6. The dotted lines in the left panels 
denote the edges deleted. In all the density plots above there is localization and therefore the EET efficiency is low. 


even a relatively small amount of disorder significantly 
increases the efficiency in an FCN. 

Finally, we have also studied the efficiency of edge- 
deleted networks in the case of N = 6 nodes with off- 
diagonal disorder and dissipation in the network, as given 
by the following term in the master equation: 

N 

LdissP = E^7n[2o-“pcrt - (H) 

n=l 

We have considered a large number of networks (28 in 



t t 


FIG. 5: Sink population versus time (in logarithmic scale) in 
the absence (solid line) and presence of dephasing with differ¬ 
ent values of the dephasing coefficient 7deph = 0.01 (dashed- 
dotted line), 0.1 (dashed line) and 1.0 (dotted line), for (a) 
a fully connected network, and networks with deletion of the 
edges (b) 1 — 6, (c) 2 — 4 and (d) 3 — 6. It can be seen that 
in the cases with localization, the influence of noise is con¬ 
structive and reduces the saturation time, while for the case 
without localization (i.e. panel (b)) noise increases the satu¬ 
ration time. Thus noise has a constructive influence on the 
cases with localization. 



FIG. 6: (Golor online) Dependence of the saturation time Ts 
on the hopping rate Jie. The saturation time is minimized 
by Jre ~ 3.04. 


total), ranging from the case where there are six con¬ 
nections between the nodes up to the FCN, with fixed 
7n = 0.01 for all sites, as a function of the disorder 
strength x between 0 < x < 0.3. While the detailed 
dependence of r^oo on x is complicated and depends on 
the topology of the network, we find that for networks for 
which rjoo (x = 0) is close to unity to start with, increasing 
X does not have much effect, as expected. In some cases, 
however, there’s a slight approximately linear decrease in 
the efficiency, up to about 10 % from x = 0 to x = 0.3. 
On the other hand, for networks where r]oo{x = 0) is low 
(less than 0.4), there’s an approximately linear increase 
in the efficiency up to about twice its value at zero disor¬ 
der (e.g. in the case of the FCN). More detailed results 
of these studies will be published elsewhere. 
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VIII. SUMMARY AND CONCLUSIONS 

In summary, we have demonstrated in this work that 
quantum transport of an energy excitation in a sym¬ 
metric network of fully connected two-state objects, de¬ 
scribed by a time reversal symmetric Hamiltonian, is 
highly inefficient. Whenever an excitation is injected into 
a node which is not directly connected to the sink, a 
wave packet travelling through closed loop paths arrives 
in-phase to the initial position due to the time reversal 
symmetry, which leads to destructive interference of the 
transition amplitudes at the initial site. The correspond¬ 
ing localization of the excitation energy wave packets in¬ 
side the network is the reason for inefficiency of transport 
in a symmetric, fully connected network. Reducing the 
symmetry by eliminating hopping between nodes not di¬ 
rectly connected to the sink preserves the localized states 
and hence does not significantly increase the system ef¬ 
ficiency. We show that an efficient way to improve the 
efficiency of energy transfer is to introduce asymmetry 
between the initial node and the one directly linked to 
the sink. The rate of energy transfer can be minimized 
by a nontrivial value of the hopping rate Jiq. We have 
also included dephasing to the network and shown, as ex¬ 
pected, that the presence of a dephasing noise increases 
the system efficiency. The same effect can be obtained by 
introducing off-diagonal disorder in the hopping integral 
parameters between the nodes. Finally, we note that we 


have also carried out similar calculations for the cases 
N = b and V = 15 to verify the conclusions presented 
here. Our results should be useful in designing quantum 
networks for novel applications. 

Acknowledgements: T.A-N. has been supported in 
part by the Academy of Finland through its COMP CoE 
grant no. 251748. We wish to thank Leonardo Novo for 
useful comments. 



X 


FIG. 7: System efficiency versus disorder strength 


^ B. Gemeinholzer, Phylogenetic Networks, in Analysis of 
Biological Networks (eds B.H. Junker and F. Schreiber), 
John Wiley & Sons, Inc., Hoboken NJ, USA, (2008). 

^ R. Pastor-Satorras, M. Rubi, and A. Diaz-Guilera, Statis¬ 
tical mechanics of complex networks (Vol. 625). Springer 
Science & Business Media (2003). 

^ Z. Zimboras, M. Faccin, Z. Kadar, J.D. Whitfield, B.P. 

Lanyon, and J. Biamonte, Sci. Rep. 3, (2013). 

^ F. Garuso, New J. Phys. 16, 055015 (2014). 

^ S. Perseguers, M. Lewenstein, A. Acm, and J.I. Girac, Nat. 
Phys. 6, 539 (2010). 

® B. Giese, J. Amaudrut, A.K.Kohler, M. Spormann, and S. 

Wessely, Nature 412, 318 (2001). 

^ H. Lee, Y.G. Gheng, and G.R. Fleming, Science 316, 1462 
(2007). 

® T. Brixner, J. Stenger, H.M. Vaswani, M. Gho, R. E. 

Blankenship, and G.R. Fleming, Nature 434, 625 (2005). 

^ G.S. Engel, T.R. Galhoun, E.L. Read, T. K. Ahn, T. Man- 
cal, Y.G. Gheng, R.E. Blankenship, and G.R. Fleming, Na¬ 
ture 446, 782 (2007). 

A. Ishizaki and G.R. Fleming, Proc. Natl. Acad. Sci. 106, 
17255 (2009). 

R.E. Fenna and B.W. Matthews, Nature 258, 573 (1975); 
J. Adolphs and T. Renger, Biophys. J. 91, 2778 (2006). 

G. Panitchayangkoon, D. Hayes, K.A. Fransted, J.R. 
Garam, E. Harel, J. Wen, R.E. Blankenship, and G.S. En¬ 
gel, Proc. Natl. Acad. Sci. USA 107, 12766 (2010). 

E. Gollini, G.Y. Wong, K.E. Wilk, P.M.G. Gurmi, P. 


Brumer, and G.D. Scholes, Nature 463, 644 (2010). 

R. Hildner, D. Brinks, J.B. Nieder, R.J. Gogdell, and N.F. 
van Hulst, Science 340, 1448 (2013). 

A. Olaya-Gastro, G.F. Lee, F.F. Olsen, and N.F. Johnson, 
Phys. Rev. B 78, 085115 (2008). 

F. Garuso, A.W. Ghin, A. Datta, S.F. Huelga, and M.B. 
Plenio, J. Ghem. Phys. 131, 105106 (2009). 

A.W. Ghin, A. Datta, F. Garuso, S.F. Huelga, and M.B. 
Plenio, New J. Phys. 12, 065002 (2010). 

M. Mohseni, P. Rebentrost, S. Lloyd, and A. Aspuru- 
Guzik, J. Ghem. Phys. 129, 174106, (2008). 

P. Rebentrost, M. Mohseni, 1. Kassal, S. Lloyd, and A. 
Aspuru-Guzik, New J. Phys. 11, 033003, (2009). 

Y. Aharonov, L. Davidovich and N. Zagury, Phys. Rev. A 
48, 1687 (1993). 

O. Miilken and A. Blumen, Phys. Rep. 502, 37 (2011). 

G. H. Bennett and G. Brassard, Proceedings of IEEE In¬ 
ternational Gonference on Gomputers, Systems and Signal 
Processing, pp. 175-179 (1984). 

G. H. Bennett, G. Brassard, G. Grepeau, R. Jozsa, A. Peres, 
and W.K. Wootters, Phys. Rev. Lett. 70, 1895 (1993). 

A. Acm, J.I. Girac, and M.M. Lewenstein, Nat. Phys. 3, 
256 (2007). 

S. Perseguers, J.I. Girac, A. Acm, M. Lewenstein, and J. 
Wehr, Phys. Rev. A 77, 022308 (2008). 

H. J. Kimble, Nature 453, 1023 (2008). 

N. Gisin and R. Thew, Nature Photonics, 1, 165 (2007). 
D.I. Tsomokos, M.B. Plenio, 1. de Vega and S.F. Huelga, 









7 


Phys. Rev. A 78, 062310 (2008). 

S. Bose, A. Casaccino, S. Mancini, S. Severini, Interna¬ 
tional Journal of Quantum Information 7, 713 (2009). 

L. Novo, S. Chakraborty, M. Mohseni, H. Neven, Y. Omar, 
Sci. Rep. 5, 13304 (2015). 

A. Barrat, M. Berthelemy, and A. Vespignani, Dynami¬ 
cal processes on complex networks, Cambridge University 
Press (2008). 

M. B. Plenio and S.F. Huelga, New J. Phys. 10, 113019 
(2008). 

H. Bassereh, V. Salari, and F. Shahbazi, J. Phys.: Cond. 


Mat. 27, 275102 (2015). 

A. Asadian, M. Tiersch, G.G. Guerreschi, J. Cai, S. 
Popescu, and H.J. Briegel, New J. Phys. 12 , 075019 (2010). 
J.R. Johansson, P.D. Nation, and F. Nori, Comp. Phys. 
Comm. 184, 1234 (2013). 

Roman S. Ingarden, A. Kossakowski, and M. Ohya, Infor¬ 
mation Dynamics and Open Systems: Classical and Quan¬ 
tum Approach, Springer Science & Business Media (2013). 
A. Rivas and S.F. Huelga, Open Quantum Systems. Berlin: 
Springer (2012). 



